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Abstract — The  rapid  growth  in  wireless  communication 
technology  has  led  to  a scarcity  of  spectrum.  But,  studies 
are  saying  that  licensed  spectrum  is  underutilized. 
Cognitive  Radio  Networks  ( CRNs)  seem  to  be  a promising 
solution  to  this  problem  by  allowing  unlicensed  users  to 
access  the  unused  spectrum  opportunistically.  In  this 
paper  we  proposed  a novel  spectrum  sensing  algorithm  to 
improve  the  probabilities  of  detection  and  false  alarm  in 
a CRN,  using  the  traditional  techniques  of  energy  and 
first  order  correlation  detection.  Results  show  a 
significant  improvement  in  performance  in  cooperative 
spectrum  sensing. 

Keywords — cognitive  radio  networks,  energy  detection, 
correlation  detection,  cooperative  spectrum  sensing. 

I.  INTRODUCTION 

CRNs  have  been  emerging  as  a promising  solution  for  the 
shortage  of  spectrum  in  wireless  communication.  The 
design  of  any  CRN  is  based  on  opportunistic  usage  of 
licensed  spectrum  without  causing  interference  to  the 
primary  user  (PU)  or  licensed  user.  Spectrum  sensing  is 
the  first  step  of  the  working  of  a cognitive  radio  network. 
The  others  being  analysis,  decision  and  adapting. 

Spectrum  sensing  is  broadly  classified  as  three  types; 
non-cooperative,  cooperative  and  interference  based.  We 
use  the  non-cooperative  techniques  of  energy  and 
correlation  detection  techniques  to  propose  an  algorithm 
for  cooperative  spectrum  sensing.  Both  energy  and 
correlation  detection  uses  a threshold  to  make  a decision 
on  whether  the  spectrum  is  occupied  byPU  or  free  for  the 
use  of  secondary  user  (SU)  or  cognitive  radio  user. 

A traditional  energy  detection  is  proposed  in  [2]  and  the 
assumption  in  the  technique  is  that  the  energy  of  the 
signal  is  higher  than  the  threshold,  if  PU  is  present.  This 
information  is  used  to  make  a decision  by  the  SU. 
Authors  in  [5]  used  covariance  information  of  received 
samples  to  determine  the  presence  or  absence  of  PU.  The 
threshold  is  generated  by  using  a covariance  matrix  of 
received  samples  by  assuming  the  desired  sample  gives  a 
higher  covariance  if  PU  is  present.  This  threshold  is  used 
to  make  a decision  for  the  future  received  signals. 

In  this  paper  we  propose  a solution  for  cooperative 
spectrum  sensing  using  correlation  based  detection. 
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Statistics  we  employed  to  measure  the  improvement  are 
probability  of  detection  (Pd)  and  probability  of  false  alarm 
(Pf).  Probability  of  detection  being  the  chance  of  SU 
detecting  the  presence  or  absence  of  PU  accurately  and 
probability  of  false  alarm  being  the  system  showing  that  a 
PU  is  present  even  if  the  spectrum  is  free,  rendering  the 
SU  losing  the  opportunity  to  use  the  spectrum.  We  first 
give  a novel  summary  statistic  as  the  first  order 
correlation  of  the  received  samples.  Sensing  is  formulated 
based  on  energy  and  correlation  detection.  We  derive  the 
probabilities  of  detection  and  false  alarm  for  different 
scenarios  of  spectrum  sensing  using  linear  combination  of 
the  summary  statistics.  The  linear  combination  gives  a 
lower  computational  complexity  and  closed  form 
expressions.  Also  it  can  lead  to  a new  statistic  where  the 
distance  between  the  mean  of  the  distributions  for  both 
hypotheses  tests  is  higher.  This  ensues  in  a lower  variance 
for  any  given  SNR.  The  assumption  in  this  scheme  is  that 
the  energy  and  first  order  correlation  of  the  samples  is 
known.  Because  of  this  information,  if  the  signal  samples 
are  correlated,  the  proposed  algorithm  significantly 
improves  the  performance  of  spectrum  sensing. 
Nonetheless,  if  the  signals  are  not  correlated,  then  this 
algorithm  acts  as  a traditional  energy  detection  scheme  as 
the  correlation  part  of  the  summation  becomes  zero.  The 
simulation  results  conveys  that  the  proposed  cooperative 
spectrum  sensing  algorithm  improves  both  the  probability 
of  detection  and  probability  of  false  alarm  of  a CRN. 

The  rest  of  the  paper  is  structured  as  follows.  Section  II 
we  define  the  system  model  with  energy  and  first  order 
correlation  summary  statistics.  Section  III  deals  with  the 
local  and  global  scenarios  for  spectrum  sensing  and 
derivations  of  probabilities  of  detection  and  false  alarm  in 
those  scenarios.  In  Section  IV,  numerical  outcomes  are 
represented.  Conclusion  is  give  in  Section  V and  last 
section  enlists  the  references  made  to  create  this  paper. 

II.  SYSTEM  MODEL 

We  considered  a CRN  with  M users,  //(/represents  the 
absence  of  a PU  and  ///represents  the  presence  of  a PU. 
With  these  two  hypotheses  at  k,h  time  instance  where,  k 
=1,2, ...,N,  the  received  signal  xi(t)  by  the  i-th  SU  is  given 
by 
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Hq-  XL(t)  = V j(t), 

H±:  xt(t)  = vt(t)  + hiS(k),i  — 1,2,3,  (1) 

where  v,(t)  and  s(k)  denote  the  additive  white  Gaussian 
noise  (AWGN)  for  the  i-th  SU  and  the  transmitted  signal 
of  PU  at  time  k.  hi  denote  the  channel  gain  between  the 
primary  transmitter  and  the  SU  i receiver.We  assume  that 
the  exact  channel  power  gains  can  be  estimated  if  the  SUs 
know  the  primary  transmit  power.  We  assume  that  the 
transmitted  signal  s(k)  and  the  set  of  noises  {vt(k)j 
areindependent  of  each  other.  Usually  the  samples  of 
transmitted  signals  are  related  to  the  modulation  schemes 
employed,  where  adjacent  samples  of  the  transmitted 
pulses  are  correlated  with  each  other.  Thus,  the  sampling 
rate  is  chosen  in  a way  that  the  set  of  noise  samples  are 
independent  of  each  other. 

SUs  calculate  two  statistics  u,  and  r,  which  denotes 
received  energy  and  first  order  correlation  of  SU  i over  an 
interval  of  N samples,  as  follows 


w-i 


ut  = ^|Xj(/c)|2;  i = 1,2,3,.. 

. , M, 

(2) 

k= 0 

N-2 

^ xt(k  + 1 )x?(k),  i = 

1,2,3 

(3) 

fc= o 


Assuming  \xi(k)\2  for  i = 1,2,..,M  to  be  independent  and 
identically  distributed  (IID)  random  variables,  and  based 
on  the  central  limit  theorem  [6],  for  large  values  of  N,  the 
statistics  fui } are  approximately  normally  distributed  with 
means 


E[Ui\  = 


where^  = 

ui 

We  define  Es as 

N-l 


No?, 

(N  + Qo?, 

Es\hi  |2 


tfl 


(4) 


^>(fc)|2. 


(5) 


Variance  is  given  by 
(No?, 


Par[Uj]  = 


(6) 


Wo 

KN  + 2Q,  Hx 
The  same  can  be  applied  to  r,with  a slight  change  of  the 
terms  x,(k  + 1)  x * (k)  and  x,(k  + 2)  x*  (k+2)  being  not 
independent  of  each  other.  So  we  consider  them  as  two 
IID  random  variable  sets.  By  central  limit  theorem  [6],  for 
adequately  higher  values  of  N,  each  of  these  sets 
approximately  follows  a normal  distribution.  As  a result, 
statistics  {rij  have  a near  normal  distribution  with  means 


E\n\  = 


o, 

\hi\2Rs(  1), 


W, 

H 


° (?) 


where 


N-2 


Rs(  1)  A y Si(k  + 1 KO), 


and  variance  is 


(8) 


Uar[rj]  = 


(9) 


(N  - 1 )o?,  H0 

[(N  -l)o?  + 2Es\hi\2o?,  H \ 

In  the  above  equation  we  make  an  assumption  that  for 
large  values  of  N,  £k=ols(^  + 1)|2  + |s(/c)|2  = 2 Es. 


The  covariance  of  M,and  r,i  s 
CovfUj.rJ  = {“]fci|2  Hh[ 


(10) 


Block  Diagram  of  System  Model 


III.  LOCAL  AND  GLOBAL  SCENARIOS  OF 
SPECTRUM  SENSING 
(a)  LOCAL  SENSING 

In  local  spectrum  sensing,  we  examine  each  SU 
individually.  Each  SU  i employs  a linear  summation  of 
u,and  r,as  a test  statistic.  It  is  defined  as 


Zi=  Wu,iUi+  wr,,r,f  1 1 ) 

where wu,i  and  w,-,  denote  the  combining  weights  for 
statistics  M,and  r„  respectively.  These  coefficients  signify 
the  contribution  of  each  statistic  to  the  decision.  Based  on 
the  above  equation  (11)  and,  equations  (4)  and  (7),  the 
mean  of  z,is 


E[zt ] = 


wU}iNo?,  H0 

wu,i(.N  + Qa?  +wri\hi\2o?Rs(l),H1 


By  (6),  (9)  and(10),  the  variance  of  zu  is  given  by 
Var[zi\Hj\  = w?EHjWi,  j = 0,1, 
wherew;=(wU  j,  wr  i)  and 


No?  0 
0 (N+l)o? 


(12) 

(13) 

(14) 


and 


zHl  - 


(15) 


(N  + 2^)o?  2\hi\2o?Rs(l) 

2\hi\2o?Rl(l)  2Es\hi\2o?  + (N-l)o? 

The  metrics  we  chose  to  measure  the  performance  of  the 
system  are  probability  of  detection  and  probability  of 
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false  alarm.  The  accurate  detection  of  the  presence  or 
absence  of  PU,  without  any  interferences,  is  detection  and 
a false  alarm  causes  a SU  to  miss  the  opportunity  to  use 
the  unused  channel.  For  local  sensing  the  decision  rule  is 
given  as 

Zi  < Yi,  H0 

Zi  > Yi,  (16) 

whereyi  is  the  threshold  for  decision  for  SU  i. 

Based  on  the  above  threshold  the  probability  of  false 
alarm  according  to  [4]  is  given  by 

Pf  = P{H,\H0-)  = Q fe==\  (17) 

\JVarlZi\Ho]] 

wherePdl lUIn)  is  the  probability  that  ///is  true  i.e.,  that 
the  PU  is  active  while  in  reality  Hois  true  i.e.,  that  the  PU 
is  not  active.  The  probability  of  detection  can  be  written 


Pd=P(H1\H1)  = Q 


Yi  — E\zi\H1 


(18) 


^Varlz^HJ 

whereP(H/\Hi)  is  the  probability  that  Hi  is  true,  i.e.,  the 
PUs  are  active  and  actually  Hi  is  true,  i.e.,  the  PU  is 
active.  By  some  changes  in  equation  (17)  we  can  write  as 
Yi  = -JVar[Zi\H0]  Q~\Pf)  + E[Zi\H0\.  (19) 

By  substituting  (19)  in  (18)  we  can  write  probability  of 
detection  as 


Pd  = Q 


Wj%0w  iQ  1(Pf)-gJ  Wj 
yjVarUilH-i] 


(20) 


whereof;  = (('jCrj2,f?s(l)|/i1|2).Our  ultimate  goal  is  to 
maximize  the  probability  of  detection.  From  the  equation 
(20)  to  maximize  that,  the  value  of  /( Wi)  must  be 
minimum, 

minw.  f(wt ) (21) 

where 


f(.Wi ) = 


Wj%0W/<2  1(Pf) 


gj  Wj 


(22) 


w 


-•Hd 


It  is  to  be  noted  that  the  threshold  for  decision  directly 
relates  to  the  weighted  vector  w;for  a given  probability  of 
false  alarm  (Pf).  For  the  decision  threshold  to  be  optimum, 
equation  (21)  is  to  be  solved  over  weight  vector  w,.  The 
optimization  problem  can  be  solved  using  methods  in  [4] 


an H T71 


Etfi  - 


N diag2(o ) + diag(5u)  + 2Esdiag(o)diag(h) 
2Rs(l)diag(o)diag(h) 


of  the  CRs  in  the  network  and  decisions  are  made  by 
individual  CRs  from  that  information. 

In  this  type  of  sensing,  both  statistics  n/andr,are 
transmitted  to  the  fusion  center  trough  a noisy  control 
channel  and  a global  test  statistic  is  calculated  in  the 
fusion  center  for  the  purpose  of  detection.  The  received 
signal  at  fusion  center  can  be  written  as 

(23) 

where  u = ( ui , 112,...,  um)  and  r = (ry,  r; tm).  We  define 

nu  = (nui,  nU2,  n„3, «„m),  in  which  the  variances  are 

collected  into  the  vector  form  Su  — (8Ul,  SU2, ... , 5U  ), 

nr  = ( n,i , nr2,nr3, «,m),  in  which  the  variances  are 

collected  into  the  vector  form  8r  = (5ri,  8r, ... , 5rM).The 
received  signals  are  normal  with  means  = E\ui  \ 

and  and  E[yr.]  = E\rL\  and  variances  Far  [yu.]  = 
Far[uj]  + (i'2.  and  Far [yr.]  = F ar | rL \ + 82..  We  can 
write  the  global  test  statistic  as  a linear  combination  of  all 
the  received  signals  as 

zc  = w Tuyu  + w Tryr  = w jy,  (24) 

wherewp  (w„rwr)  in  which  w„=  (w„v  w U2,~-  w „m)  and  w,= 
(w ,1,  Wr2,~-  w,m)  are  the  weight  vectors  which  represent 
the  contribution  of  each  received  signal  in  the  decision 
making.  The  mean  of  zcis  given  by 


E[zc ] = 


Nw  u<7, 


(25) 


(w u(No  + Es h)  + wr  hRs(l),  H f 
where  h =(\hi\2,\h2\2,...MiM\2),  o =(oi2,022,...,Om2). 

The  variances  under  different  hypotheses  are  given  by 
Var[zc\Hi\  = wjz„.wf,  i = 0,1,  (26) 

whereE/)oand  E/)i  (given  at  the  end  of  the  page  as  equation 
28)  are  given  as, 

TV  diag2(o)  + diag(8u)  [0 ]MXM 

[0]mxm  (Af  — l)diag2(o)  + diag{8r) 

(27) 

In  global  sensing  the  detector  can  be  defined  as 
Zi  < Yu  H0 

Zi  > Yi.  Hi 

Similar  to  the  probability  of  detection  of  local  sensing 


/ r 


2Rs(l)diag(o)diag  (h) 

(N  — l)diag2(o ) + diag(8r ) + 2Esdiag(o)diag(h) 


(29) 

.(28) 


CRs  are  connected  to  an  infrastructure  and,  having  the 
sensing  results  of  all  CRs,  decisions  are  made  by  the 
centralized  system  for  each  CR.  In  a de -centralized 
system,  all  CRs  are  connected  to  each  other  and  each  CR 
has  access  to  information  about  the  sensing  results  of  rest 


where^  = (Fsh, Rs(l)h).  For  a specified  probability  of 
false  alarm,  the  solution  for  maximizing  probability  of 
detection  is 


mm 


Wf 


f(Wf) 


(31) 
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where 

w/XHow/  Q-1^)  - gTwf 

fW  = , (32) 

Jw  JlHlwf 

This  optimization  problem  can  also  be  solved  by  the 
proposed  methods  in  [4]  and  [7]. 

IV.  RESULTS  AND  GRAPHS 

In  this  section  we  simulate  the  scheme  for  different 
scenarios.  We  assumed  M CR  users  in  the  network 
topology.  For  reducing  the  complexity  we  have  taken  the 
transmitted  primary  signal  as  s(k)  =1.  In  Fig.  1,  for  the 
simulation,  we  considered  a cognitive  radio  network  with 
M = 6,  N = 20,  or  = 1,  Vi.  i = l,Vi.  We  took  the  local 
SNRs  at  each  CRs  as  [9.3,7.8,9.6,7.6,3.5,9.2]  in  dB.  The 
results  we  achieved  are  based  on  106  noise  realizations. 
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Fig.l : Single  CR  v.v  Multiple  CRs  for  the  proposed  model 

For  Fig.  2,  the  simulation  parameters  are  N = 20,  5Ui  = 6,-i 
= l,Vi  and  the  graph  is  depicted  for  values  of  M = 1,3,6 
and  10,  where  the  corresponding  values  of  average  local 
SNRs  are  9. 3,6. 7, 9. 2 and  9.0  in  dB.  For  M = 1,  the  noise 
level  is  considered  as  a = 1.9  and  the  local  SNR  is  8.3. 
For  M = 3,  the  noise  level  is  c = {0.7,  1.0,  0.8}  and  the 
local  SNRs  are  {10.4,  9.3,  2.6}dB.  For  M = 6,  the  noise 
levels  are  considered  as  o={0.9,  1.3,  1.0,  2.0,  1.8,  1.2} 
and  local  SNRs  are  {7.2,  5.1,  10.8,  -1.2,  3.6,  9.7}  in  dB. 
For  M = 10,  thenoise  levels  are  considered  as  o={0.9,  1.3, 
1.0,  2.0,  1.8,  1.2,  1.5,  0.8,  1.1,  1.8}  and  local  SNRs  are 
{7.2,  5.1,  10.8,  -1.2,  3.6,  9.7,  9.0,  8.4,  6.0,  -0.8}  in  dB. 
The  probabilities  in  Fig.2  clearly  shows  that,  as  the 
number  of  CRs  in  the  network  increases,  the  probability 
of  detection  increases  and  the  probability  of  false  alarm 
decreases  because  of  the  cooperation  between  the  CRs. 
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Fig.2:  Probability  of  detection  vs  probability  of  false 
alarm  under  the  proposed  approach  for  number  of  users 
M =1.3,6,10. 

V.  CONCLUSION 

In  this  paper  we  have  proposed  a novel  algorithm  for 
cooperative  spectrum  sensing  in  cognitive  radio  networks 
to  increase  probability  of  detection  and  decrease 
probability  of  false  alarm.  We  have  used  two  methods  of 
non-cooperative  spectrum  sensing,  energy  and  first  order 
correlation  detection  and,  the  algorithm  has  been 
formulated  by  the  linear  combination  of  both.  The 
simulations  show  that  the  algorithm  improves  by  number 
users  in  the  network  and  is  an  improvement  on  the 
traditional  non-cooperative  methods. 
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